Calculation of conduction-to-conduction and valence-to-valence transitions between 
bound states in (In,Ga)As/GaAs quantum dots 
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We have calculated the conduction-to-conduction and valence-to- valence absorption spectrum of 
bound states in (In,Ga)As/GaAs quantum dots charged with up to three electrons or holes. Sev- 
eral features emerge: (i) In pure (non-alloyed) InAs/GaAs dots, the lS-lPi and IS-IP2 conduction 
intraband transitions are fully in-plane polarized along [110] and [110], respectively, while valence 
transitions are weakly polarized because the hole P states do not show any in-plane preferential 
orientation, (ii) In alloyed Ino.6Gao.4As/GaAs dots the [110] and [110] polarization of the corre- 
sponding 1S-1P conduction intraband transitions is weakened since the two IP states are mixed by 
alloy fluctuations. The polarization of valence intraband transitions is insensitive to changes in alloy 
fluctuations, (iii) For light polarized along [001], we find a strong valence-to- valence transition that 
involves a weakly confined hole state with predominant light-hole character, (iv) When charging 
the dots with a few electrons, the conduction intraband transitions display spectroscopic shifts of 
~ 1-2 meV. These shifts are a result of correlation effects (captured by configuration-interaction) 
and not well described within the Hartree-Fock approximation, (v) When charging the dots with 
holes, valence intraband spectra are more complex than the conduction intraband spectra as hole 
states are strongly affected by spin-orbit coupling, and configuration mixing is more pronounced. 
Spectroscopic shifts can no longer be identified unambiguously. These predictions could be tested 
in single-dot spectroscopy of n-doped and p-doped quantum dots. 



I. BRIEF BACKGROUND AND SCOPE 

Doping quantum dots w-type (p-type) followed 
by infrared-light excitation leads to conduction-to- 
conduction (valence-to- valence) intraband excitations be- 
tween confined states. In early experiments, intraband 
transitions were studied for quantum wells^ where in 
the case of n-type doping the selection rules lead to 
allowed transitions only for light polarized normal to 
the well^i 3 -^ whereas p-type samples afforded transi- 
tions allowed in normal incidence. Instead, 77,-doped 
quantum dots offered the possibility of allowed intra- 
band transitions for light incident normal to the plane 
of the dots. Thus, quantum dots became the focus 
of much experimental work devoted to study the in- 
traband optical transitions. 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ! 12 ^ 3 ! 14 ! 15 For ex- 
ample, measurements of optical intraband transitions in 
(In,Ga)As/GaAs quantum dots charged with electron 
and holes have been recently performed by Zibik et al*£, 
who measured conduction-to-conduction transitions, and 
Preisler et a/pi 7 - who have measured valence-to- valence in- 
traband absorption. 

On the other hand, calculations of the intraband op- 
tical properties are often based on mode l 18 ! 19 ' 20 or k • 
p 21 approaches and not high-level atomistic approaches. 
Here, we address this issue by calculating conduction- 
to-conduction and valence-to-valence intraband optical 
absorption of (In,Ga)As/GaAs dots charged with up to 
three electrons or holes. We do not intend to survey 
the effects of size, shape, and composition on the in- 
traband transitions of the dots. We focus primarily on 
(i) the effects of alloy fluctuations on the polarization 
of the intraband transitions, and (ii) the differences be- 
tween the absorption spectra of n-doped and p-doped 



dots as well as the spectroscopic shifts induced by charg- 
ing. We use a combined approach to calculate the in- 
traband absorption spectra in which we find the electron 
and hole single-particle states of the dots with an atom- 
istic pseudopotential method and solve the many-particle 
states of charged dots within a configuration-interaction 
approach. The advantage of this approach over simpli- 
fied methods is that it naturally includes (i) the correct 
symmetry of the dot; (ii) strain and alloying effects; (iii) 
multiband and multivalley coupling; and (iv) spin-orbit 
interaction. Our approach is purely electronic and does 
not include polaron (electron-phonon coupling) effects. 
We illustrate our findings with a prototypical lens-shaped 
Ino.6Gao.4As/GaAs dot with diameter b — 252 A and 
height h = 35 A. As a benchmark, for dots charged with 
a single carrier, we provide results for a pure non-alloyed 
InAs/GaAs dot with the same size. While so far in- 
traband experiments have focused on ensemble of dots 
with different degrees of homogeneity, which broadens 
the observed transitions, we offer predictions that could 
be probed in single-dot spectroscopy of n- or p-doped 
dots. In addition, as the control of doping carriers is 
difficult, i.e in an ensemble of doped dots some have a 
single electron while other have two or none, we discuss 
the intraband spectra upon increasing the carriers one at 
a time. 



II. METHOD: PSEUPOTENTIAL APPROACH 
AND CONFIGURATION-INTERACTION 

We calculate the single-particle electron and hole en- 
ergy levels of the self-assembled dot within an atom- 
istic, pseudopotential-based method^ The wave func- 
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tion ipj and energy Ej are solutions of the atomistic 
single-particle Schrodinger equation 



Vso+Y. «a(R~Rl,a) 



i']= £ jip3 (!) 



where v a is the screened pseudopotential for atom of type 
a (In,Ga,As) with position R; a within the dot or barrier, 
and Vso is a non-local pseudopotential that accounts for 
the spin-orbit interaction^ These pseudopotentials are 
carefully fitted to bulk GaAs, InAs, and (In,Ga)As alloys, 
thus removing the local-density-approximation (LDA) er- 
rors. The basis in which we expand ipj to solve Eq. fl} is 
a linear combination of full-zone Bloch bands of the un- 
derlying solids, 24 Thus, this method incorporates multi- 
band and multivalley coupling, band non-parabolicity, 
and spin-orbit effects, as well as the effects of the un- 
derlying strain in the dot and barrier. 

To solve for the many-particle states {^f v (Af), E U (N)} 
of the dot with M carriers, where Af = N e electrons 
or Nh holes, we use a configuration-interaction approach 
with screened direct (J) electron-electron and hole-hole 
Coulomb interaction and exchange (K)^ This method 
has been recently applied to the calculation of electronic 
and optical properties of (In,Ga)As/GaAs dots such as 
electron and hole charging 2 ^, radiative lifetimes of neu- 
tral and charged excitons^ relaxation times of elec- 
trons due to electron-hole Auger scattering,- 29 and fine- 
structure splittings of neutral and charged excitons^ 

At low temperatures such that only the ground state 
i S?o(AT) of the Af-carrier dot is significantly occupied, the 
optical absorption spectrum for light polarized along e is 
given by 



I(hu,;e) = J2\( 1 *»'\e-P\yo)\ 2 5(E lJ , - E -huj). (2) 



In the results we present subsequently, we have phe- 
nomenologically broadened the spectra with a Gaussian 
of width a — 0.25 meV. Such a broadening has also 
been used in other simulations of optical absorption in 
the literature. 

For dots with cylindrical symmetry, in which the pro- 
jection of the angular momentum along the cylindrical 
axis (z) is a good quantum number, the selection rules 
for intraband transitions are the following. For e || [100] 
or [010] transitions are allowed between states such that 
AL Z = ±1, while for e || [001] only transitions with 
AL Z = are allowed. 

We consider in our calculations prototypical lens- 
shaped Ino.6Gao.4As/GaAs dots with base diameter b = 
252 A and height h = 35 A. As a benchmark, in the 
case of dots charged with a single carrier we also con- 
sider a pure non-alloyed InAs/GaAs dot of the same 
sizej^l Note that the detailed experimental character- 
ization of the shape, size, and composition of the al- 
loyed (In,Ga)As/GaAs dots probed optically is scarce. 



Electron and hole energy levels 
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FIG. 1: Electron (top) and hole (bottom) energy levels in lens- 
shaped (In,Ga)As/GaAs dots with base diameter b = 252 A 
and height h = 35 A relative to the GaAs barrier conduc- 
tion band minimum _E c (GaAs) = —4.093 eV (w.r.t vacuum) 
and its valence band maximum _B„(GaAs) = —5.620 eV 
(w.r.t vacuum). The energy levels are more confined in the 
pure (non-alloyed) InAs/GaAs dot (a) than in the alloyed 
Ino.6Gao.4As/GaAs dot (b). For holes, we show the first 40 
confined states. Black and grey arrows indicate, respectively, 
the lowest intraband conduction transition for polarization 
along [110] and [110]; while numbers show the transition en- 
ergy in meV. 



The prototypical alloyed dot we consider here presents 
properties like excitonic gap; electron and hole single- 
particle energy spacings; and binding energies of neutral 
and excited excitons in excellent agreement with avail- 
able data^ 



III. RESULTS FOR SINGLE-PARTICLE 
ELECTRON AND HOLE LEVELS 

A. Electron levels 

Bound states can be labeled by their leading or- 
bital character and approximately arranged into shells: 
{IS; lPi, 1P 2 ; 1-Di, 1D 2 , 25; . . .}. For the dot size con- 
sidered in our calculation, a pure (non-alloyed) dot con- 
fines 15 states arranged in five shells [Fig. QJa)], two 
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more shells than in its alloyed counterpart which binds 
10 states [Fig. (Hb)]. The larger number of confined lev- 
els in a pure dot is due to a larger strain-modified con- 
duction band offset than that in an alloyed dot. In the 
pure InAs/GaAs dot [Fig. QJa)] the 15 state is located 
306 meV below the conduction band minimum (CBM) of 
the GaAs barrier, while in the alloyed dots [Fig. HJb)] is 
-118 meV below the GaAs CBM. (This value changes 
by a few meVs depending on the random alloy fluctua- 
tions in the dot.) These energies of the 15 state rela- 
tive to the GaAs CBM set the cutoff for conduction-to- 
conduction intraband transitions between bound states. 
The P-shell consists of two non-degenerate states lPi 
and IP2. The origin of this splitting is atomistic^ and 
a consequence of the underlying C% v symmetry of the 
(pure, nonalloyed) dots, which is lower than Coo V sym- 
metry of the macroscopic (lens) shape, which is nor- 
mally assumed in continuum effective-mass models. In 
our calculation we consider perfectly cylindrical dots; 
yet, the P-P splitting is as large as 2 meV in pure 
dots. Piezoelectricity^ and non-cylindrical shaped fur- 
ther contribute to this P-P splitting. In addition, each of 
these P states present a nearly equal mixture of L z = ±1 
components contrary to the axially-symmetric case in 
which each state has a well defined L z component. In 
the pure (non-alloyed) dot 1P\ is oriented along [110] 
and IP2 is oriented along [110] [Fig. [5(a)]. In alloyed 
dots the symmetry of the dots is lower than C^v due to 
random alloy fluctuations. In this case, not only [110] and 
[110] are mixed and different disorder realizations (fluc- 
tuations) change the magnitude of the P-P splitting by 
1-3 meV but, more remarkably, alloy fluctuations affect 
the in-plane orientation (polarization) of the P states, as 
is shown in Fig. [5J Or, equivalently, alloy fluctuations 
change the relative phase (f>± of the L z = ±1 components 
in the lPi and IP2 states, which results in different in- 
plane orientations (polarizations) of these states. For in- 
stance, in dot A we have ~ and ~ 7r/2, while for 
dot C we have 0+ ~ 0_ ~ 0. In turn, the P-shell con- 
sists of non-degenerate 1-Di, ID2, and 25 states. States 
lDi and IP2 show a nearly even mixture of L z = ±2 
components. Depending on alloy fluctuations, state 25 
can also have sizeable L z = ±2 components, thus making 
it not possible in those cases to assign a leading orbital 
character to these P-shell states. 



B. Hole levels 

Both non-alloyed and alloyed dots confine a large 
(Mh > 20) number of single-particle levels. Due to 
the multi-band nature of these hole states and for flat 
dots like the one we consider here, only low-lying states 
present shell structure that is less pronounced, i.e. larger 
P-P and D-D splittings, than in the electron case. For 
these states, one can still use their leading 5, P, D or- 
bital character to identify them. Higher lying states show 
heavy mixing of orbital character. For pure and alloyed 
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FIG. 2: (Color online) Isosurface representation of the elec- 
tron (top) and hole (bottom) wavefunctions of states lPi 
and IP2 in a pure non-alloyed InAs/GaAs dot (a), and three 
Ino.6Gao.4As/GaAs dots (b) with different alloy fluctuations. 
Dots are lens-shaped and have the same size: b = 252 A and 
h = 35 A. 



dots the 15 hole state is located, respectively, 211 meV 
and 186 meV above the valence band maximum of the 
GaAs barrier; see Fig. [TJ These values are the cutoff for 
valence-to-valence intraband transitions between bound 
states. In addition, the P states are not oriented along a 
preferential in-plane direction [Fig. [2j. 



IV. INTRA-CONDUCTION AND 
INTRA- VALENCE TRANSITIONS IN DOTS 
WITH A SINGLE CARRIER 

Assuming that only the 15 state is occupied by doping, 
for in-plane polarized light we expect intra-conduction 
transitions between bound states that satisfy AL Z = ±1 
when the dot is occupied by a single electron. The lowest- 
energy transitions correspond to 15-lPi and 15-lp? (in- 
dicated by arrows in Fig. [I]). The orientation (polariza- 
tion) of the IP states determines the polarization prop- 
erties of these transitions: 

(i) For pure (non-alloyed) InAs/GaAs dots, we expect 
that for polarization e || [110] only transition 15-lPi to 
be allowed, while for e |[ [110] only 15-1P2. 

(ii) For alloyed dots, alloy fluctuations dictates the ori- 
entation of states lPi and IP2 [Fig. |2Jb)]. Thus, we ex- 
pect a strong dot-to-dot dependence of the polarization 
properties of the two conduction intraband transitions. 
In particular, for the in-plane polarization e || [110] we 
expect both transitions 15-lPi and 15-lp2 to be allowed, 
with intensities 



4 



£ 10- 



Intraband optical absorption in alloyed dot A 
ln 06 Ga 04 As/GaAs; b=252A, h=35A 



Conduction; N p -1 



Valence; N h -1 



A, 



(a) e || [110] 



| e|| [001] | 



40 80 120 





1 10 i 

| (b) e || [1 10] | 










6 - 












4- 








- A< 


' 


2 - 
0- 




x3 


x3 

i i i iVi'i 



35 40 45 50 55 

Energy (meV) 



20 40 60 80 

Energy (meV) 



Conduction intraband optical absorption in alloyed dots B, C 
ln nE Ga„,As/GaAs; b=252A, h=35A 




35 40 45 50 55 60 35 40 45 50 55 



Energy (meV) 



Energy (meV) 



FIG. 3: (Color online) Intraband absorption spectrum for 
in-plane polarization e || [110] (a) and e || [llO] (b) in an 
Ino.6Gao.4As/GaAs dot (dot A) charged with a single elec- 
tron (left panels) and a single hole (right panels). For these 
in-plane polarizations, conduction-to-conduction transitions 
are nearly an order of magnitude more intense that valence- 
to- valence transitions (note different ordinate scales). Inset: 
Valence intraband transition for light polarized out-of-plane 
e || [001]. The strong feature corresponds to the IShh-lSih 
transition. 



(15|e • p|lPi)| 2 oc 1 - sin(0+ - (3) 
|(15|e-p|lP 2 )| 2 (x|l + sin( f > + -^)|. (4) 



A. In-plane transitions 

In a pure (non-alloyed) InAs/GaAs charged with one 
electron the lowest transition 15-1 Pi appears at 56 meV 
and is fully polarized along [110], while the second tran- 
sition 15-lp2 is located at 58 meV and is fully po- 
larized along [110]. In contrast, in alloyed dots [Figs. 
[3]Ja) and[3]^b)] we only find partial in-plane polarization: 
For dot A we find a transition at 46 and another at 
48 meV corresponding; the lowest-energy transition 15- 
lPi is partially (~ 75 %) polarized along e = [110] while 
the higher-energy 15-1P2 is partially (~ 75%) polarized 
along e = [110]. 

The valence intraband transitions have a smaller os- 
cillator strength than those of the conduction transi- 
tions. For holes, as in the case of electrons, we find two 
strong transitions corresponding to the 15^-lP^ transi- 
tions. However, these transitions are weakly in-plane po- 
larized because the hole P states do not show any pref- 
erential in-plane orientation. In addition, we also find 
weaker transitions in the interval 40-60 meV, which arise 
from the mixing between the D-shell and P-shell hole 
states. For dot A, the lowest valence intraband transition 



FIG. 4: Idem Fig. C2» and Fig. Eftb) for two 
Ino.6Gao.4As/GaAs dots with the same size as dot A but with 
different random alloy fluctuations. 

shows a higher intensity when light is polarized along the 
[110] direction, as in the experiment of Preisler el al^ 

Calculations based on simplified models are not ca- 
pable of explicitly introducing alloy fluctuations. In- 
stead, these fluctuations are naturally included within 
our atomistic approach. Figure 0] illustrates the ef- 
fect of random alloy fluctuations on the polarization 
properties of the conduction intraband transitions in 
Ino.6Gao.4As/GaAs dots with the same size as dot A but 
different random alloy fluctuations: (i) Under [110] po- 
larization, dot B presents transition 15-lPi nearly fully 
polarized and transition 15-1P2 nearly forbidden; con- 
versely, for e || [110] transition 15-lp2 is nearly fully 
polarized and 15-lPi nearly forbidden. These polariza- 
tion properties are switched when compared to dot A. 
In addition, that the lowest conduction intraband tran- 
sition in dot B is mainly polarized along [110] in both is 
in agreement with the experiment of Zibik et al^ (ii) 
In contrast, dot C presents both transitions allowed for 
polarizations [110] and [110], with very small in-plane po- 
larization anisotropics. We find that transition 15-lPi is 
polarized along [100], with transition 15-1P2 forbidden, 
and that for e || [010] transition 15-lp2 is allowed while 
15-1 Pi is forbidden. 



B. Out-of-plane polarization — intra-valence 
transitions 

The inset of Fig. [3]Ja) shows the valence intraband 
transition for e || [001]. We find a strong feature orig- 
inated from the lShh-^Sih transition, which involves a 
weakly confined, highly excited hole state with predomi- 
nant light-hole character. This transition is nearly three 
times as intense as the in-plane valence transitions. This 
transition is consistent with the selection rule AL Z = 
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Intraband optical absorption vs dot occupation 
ln 06 Ga 04 As/GaAs; b=252A, h=35A, dot B e || [1 1 0] 
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FIG. 5: (Color online) For in-plane polarization e || [110], in- 
traband conduction-to-conduction (a) and valence-to-valence 
(b) transitions calculated at the CI level versus number of 
electrons N e and holes Nh in the dot. Black, short bars in 
panel (a) for N e — 2 and 3 show the IS 1 - IP transition energies 
predicted within the Hartree-Fock approximation — Eqs. (|13p 
and (TT7]) . respectively. For N e = 2, the HF singlet (HF s) and 
the configuration-interaction singlet (CI s) and triplet (CI t) 
are indicated. The latter being forbidden. 



for this light polarization. 



V. INTRA-CONDUCTION AND 
INTRA- VALENCE TRANSITIONS IN DOTS 
WITH A FEW CARRIERS 

We now study the conduction and valence intraband 
transitions for J\F = 2, and 3 carriers occupying the dot. 

The energy of conduction and valence band transitions 
in the presence of N e electron or Nh holes is dictated by 
differences in total energies. [See Eq. Therefore, we 
expect spectroscopic shifts of the transitions upon charg- 
ing. We define the intraband spectroscopic shift Aj/ of 
transition v&j — > ^/ f upon adding a carrier as 



A. Intra-conduction transitions vs N e 

We present Hartree-Fock (HF) expressions for the in- 
traband transition energies when the dot is charged 
with N e — 1, 2, 3 and then compare to configuration- 
interaction calculations. To illustrate our findings, we 
consider dot B, which has lPi and IP2 states that are 
nearly fully polarized along [110] and [110], respectively, 
and focus on light polarized along [110]. 

N e — 1: The ground state is |^o) = l e o) an d under 
IR light the final state is \e\). (Both states are twofold 
degenerate.) For this occupation, there are no many- 
body effects and the energy of the initial and final state 
are simply given by 



£0(1) - £} e) 



s 



E f {i) = 

The intraband transition energy is 
hw if (l)=E f -E i = 



(7) 
(8) 



(9) 



N e = 2: The closed-shell (nondegenerate) state |^o) = 
ep) is the ground state, and there are four possible final 
states originating from leje}). At the single-particle level 
the four final states are degenerate, but within the HF 
approximation these states split in a triplet (t) and a 
singlet (s): |eje})t and |ejej) s . The splitting between 
these states is due to the electron-electron exchange in- 
teraction, between an electron in \S) and an electron in 

I Pi), and equals 2Kgp. Under IR light excitation, only 
the singlet is an allowed final state. Thus, we have the 
following energies for the initial and final state: 



Eq(2) = 2£$+J t 



.(ee) 

ss 



Ef{2) 



°is 



J. 



S (2) 

(ee) 
SPi 



K 



(ee) 
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(10) 
(11) 



This results in a IS- IP transition energy in the presence 
of an extra electron that can be written as 



Ai,(.A0 = fo*/(A0 - hu>if(J\f - 1). (5) 

The energy of a many-particle state can be ex- 

pressed as superposition of a Hartree-Fock term [E ^ HF ^] 
and correlation component (5); namely, 

Ei CI \M) - El HF \N) + <y„(A0. (6) 

For the conduction-to-conduction intraband transition 
energies and shifts, we compare results for the conduction 
intraband transition energies and shifts obtained within 
the Hartree-Fock approximation and CI calculations. 



7^(2) = £l;l-£^+J^+K^^4p + [5 f (2y5 (2)]. 

(12) 

Here, Jg^ and J^spI are J respectively, the direct 
Coulomb interaction between two electrons in IS), and 
one electron in \S) and another in |Pi). For dot B our 
atomistic calculation predicts Jgg^ — 20.8 meV and 

j sp! = 16 - 6 mcV > and K SPi = 4 -3 meV. 

Our CI calculations for N e — 2 are shown in Fig. OHa). 
We find that the four states arising from configuration 
eje} indeed split in a singlet and triplet and, as expected, 
only the intraband transition from the ground-state to 
the singlet is allowed — strong feature at ~ 45 meV. The 
Hartree-Fock predicted transition energy [first five terms 



6 



in Eq. (fl2|) ] for dot B appears as a short black bar in Fig. 
[5]Ja) and is 43.7 meV, while CI predicts 44.9 meV, so the 
correlation correction [terms in brackets in Eq. (|12p ] is 

1.2 meV. 

The spectroscopic shift between the transitions in 
N e = 1 and 2 is thus given by 

A(2) = J$-4?+K$ + [8 f (2)-8o(2)](13) 

The Hartree-Fock component of this shift is 0.1 meV, 
whereas our CI calculation reveals a blue-shift of A (2) = 

1.3 meV; i.e., 1.2 meV from correlation effects [terms 
in bracket in Eq. (fT3"|) ]. This shows correlation effects 
determine the spectroscopic shift. 

N e — 3: The ground state of the initial state is |e§e}) 
(twofold degenerate) and in addition to the final states 
originating from the P shell: \e\e\) and leje^), one now 
has six states arising from transitions to the D shell: 
e o e 3)> l e o e i)> an d |e§eg) . The energy of the ground state 
and the P-shell derived final are the following: 

E {3) = 2£^ +£$+4? +24P -K^f + <^I) 
E fP (3) = ie[% + £$ +2 4 e p e) +J P p + S f (3). (15) 

Here, J^p is the direct Coulomb interaction of two elec- 
I 



We find in our CI calculations [Fig. EJa)] that the strong 
feature around 49 meV corresponds to the (nearly over- 
lapping) IP-lDx and IP-ID2 transitions. The weak 
transition at ~ 52 meV arises from 1P1-2S* and because 
\2S) in the alloyed dot is primarily oriented (polarized) 
along [110] the transition is weak. 



B. Valence transitions vs Nh 

Earlier calculations assumed simple models with the 
incorrect symmetry and neglected the multiband nature 
of the hole single-particle states, which leads to an incor- 
rect treatment of the hole-hole interaction. Within our 
atomistic approach, spin-orbit coupling and the multi- 
band nature of the hole single-particle states prevent us 



trons occupying \Pi). For dot B, J P p = 16.0 meV. The 
resulting 1S-1P\ intraband transition energy is 



7*0,0/(3) = £t P \-4^+4{%-4f+K^ i+ [5 f {i)-5 {Z)]. 

(16) 

For dot B, the HF component amounts to 43.1 meV 
[black bar in Fig. [5ja)] and the term due to correlations 
(in brackets) amounts to 2.9 meV. 

The spectroscopic shift between the lS'-lPi transitions 
in N e = 2 and 3 is 



A(3) = J P e X - 4p! + [*/(3) - S f (2) - <5 (3) + *o(2)] . 

(17) 

The HF part is a red-shift of A(3) = -0.6 meV. The 
value predicted by EMA is nearly twice as big; namely, 
[A(3)]ema = —1-3 meV. Yet, our CI calculation shows 
[Fig. \5{a) for N e = 3] that the transition is actually 
blue-shifted by 1.1 meV with respect to that in N e — 2. 
The correlations contribution [in brackets in Eq. (|17|) ] 
being 1.7 meV. Note also that the intensity of transition 
15-lPi is bleached (reduced) by nearly 50% as a conse- 
quence of Pauli blocking — one electron occupies |1Pi). 

Regarding the P-D transitions, within the HF approxi- 
mation we expect them to occur at the following energies: 



(18) 
(19) 
(20) 



from writing meaningful HF expressions in the case of 
holes. So we discuss directly the results of our CI cal- 
culations. Figure [5^b) shows the valence intraband tran- 
sitions for Nh — 1, 2, and 3 for light polarized along 
e || [110]. In general, compared to the conduction case, 
the valence intraband spectra are more sensitive to the 
number of holes in the dot. 

Nh = 2: The lS-lPi transition (lowest feature in 
Nh = 1) appears red-shifted by nearly 6 meV and split — 
two peaks between 8-11 meV. Due to the hole-hole ex- 
change interaction this transition is split in a pair of low- 
energy, nearly doubly degenerate states and two higher- 
lying states mutually split by ~ 1 meV. Similarly, tran- 
sition 1S'-1P2 splits in two transitions: One transition 
at ~ 14 meV, with an ensuing red-shift of 6 meV, and 
another at ~ 20 meV that appears slightly blue-shifted 



fc*(3) = £[% - £[% + 2 - K&l - 2 4%\ - 4g + 5,(3) 

^ 2 ( 3 ) = e[% - s[% + 2 4-) - 4 e ;> - 2 43 - 45 + ^(3) 

7k,3(3) = - 4% + 2 4") - 4") - 2 - 4 e 2 e ] + S 3 (3) 
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(~ 1 meV) with respect to the transition at Nh = 1. 

Note that contrary to the case of electrons, and due 
to the presence of spin-orbit interaction, the four states 
arising from the two- hole configuration h^hl, do not split 
in a triplet and one singlet. Instead, these four states 
split in two doublets that are allowed under IR light ex- 
citation. More importantly, in the commonly used EMA 
with two-dimensional harmonic confinement and without 
spin-orbit coupling one would not find these double-peak 
structure of allowed transitions, but instead one would 
find a spectra that resembles that of the N e — 2 electron 
case. 

Nh = 3: While the ground state is well described 
by |/ig/i}}, the effect of configuration mixing in the final 
states (upon absorption) due to hole-hole interaction is 
remarkably pronounced and leads to a complex spectrum. 
As a result, it is not possible to determine unambiguously 
the spectroscopic shifts A(3). Prominent features are the 
following. 

(i) The lowest-energy peak at nearly 9 meV corre- 
sponds to transition IS-lPi. Also, the peak at ~ 15 meV 
is mainly IS-lPi, but mixed with 1Pi-2S and lPi-25. 
This mixing leads to the high intensity of this transition. 
Remarkably, we find that in contrast to the N e = 3 case, 
the IS-lPi transition is not bleached by having a hole 
occupying the 1P\ state. 

(ii) The peak at 10 meV correspond to transition IS- 
IP2, while the weaker feature at ~ 11 meV is due to a 
transition with IP1-ID1 predominant character. 

(iii) Above 15 meV the features in the spectrum corre- 
spond to transitions to heavily mixed final configurations: 
(a) The lower-energy peak in the double-peak structure 
around 20 meV corresponds to a mixture of the allowed 
lS'-lPi and IS-IP2 transitions in addition to a sizeable 
component (16%) of the forbidden hghl-h^h^ transition. 
In turn, the higher-energy peak is a mixture of IS-IP2 
and \P\-Di transitions, (b) The peak at 26 meV arises 
from two nearly overlapping transitions. These transi- 
tions are a mixture of allowed P-D and P-F transitions, 
as well as forbidden transitions. 

(iv) Although significantly weaker than the other fea- 
tures in the spectrum, the peak at 24 meV corresponds 
to a forbidden transition made allowed by configura- 
tion mixing with allowed transitions. We also have 
found this type of transitions in the interband spectra 
of (In,Ga)As/GaAs dots^ 



VI. SUMMARY 

By combining an atomistic, pseudopotential-based ap- 
proach with the configuration method, we have cal- 
culated the conduction and intraband transitions in 
(In,Ga)As/GaAs quantum dots with up to three carriers. 
We illustrated our calculations with a prototypical lens- 
shaped Ino.6Gao.4As/GaAs dot with diameter b = 252 A 
and height h = 35 A. And as a benchmark, for dots 
charged with a single carrier, we provided results for a 
pure non-alloyed InAs/GaAs dot with the same size. We 
have made specific predictions that could be probed in 
single-dot infrared spectroscopy of n-doped and p-doped 
dot: 

(i) In pure, non-alloyed InAs/GaAs dots, the 1S-1P 
conduction intraband transitions are fully in-plane po- 
larized, while valence transitions are weakly polarized 
because the hole P states do not show any in-plane pref- 
erential orientation. 

(ii) In alloyed Ino.6Gao.4As/GaAs dots the in-plane 
polarization of IS-1P conduction intraband transitions 
strongly depend on alloy fluctuations, which change the 
in-plane orientation of the nearly generate P-shell states. 
The polarization of valence intraband transitions is in- 
sensitive to changes in alloy fluctuations. 

(iii) Upon changing the number of carriers in the dot, 
the intraband transitions display spectroscopic shifts of 
about 1-2 meV. These shifts are not well described 
within Hartree-Fock, instead their magnitude is deter- 
mined by correlation effects. 

(iv) Spin-orbit coupling and the multi-band character- 
istic of holes states result in important differences be- 
tween the many-particle valence and conduction intra- 
band spectra. Spectroscopic shifts can only be deter- 
mined unambiguously for conduction transitions. 
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